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: 1 Introduction 

J3 • The Kepler problem is an ancient problem, solved for the first time by Newton 

and well studied over more than three centuries jTTJ, [9], [7], [3], [2]. Newton 
showed that the solutions of his equation of motion in Hamiltonian form 



q = p , p = -q/<? 3 

are planar ellipses traversed according to the area law. What is commonly 
called the Kepler problem is a study of the solutions of these equations, 
preferably from a geometric point of view and taking symmetry arguments 
into account. It might come as a surprise that the past half century we have 
q ■ still witnessed new insights on the Kepler problem. 

A key point in a modern treatment of the Kepler problem in M. n is the 
regularization of the collision orbits, which goes back to Moser [10]. The 
treatment of Moser relates the Kepler flow for a fixed negative energy level 
to the geodesic flow on the sphere § n . The method of Moser is geometrically 
well motivated. The main calculation to be done is the determination of the 
explicit transformation formulas for stereographic projection on the level of 
cotangent bundles. 

An alternative approach to the regularization, having the advantage of 
transforming canonically the negative energy part of phase space R 2n for the 
Kepler Hamiltonian into the punctured cotangent bundle of S n , was found by 
Ligon and Schaaf [8] . Unfortunately the article by Ligon and Schaaf requires 
a good deal of computations. Cushman and Duistermaat have given a simpli- 
fied treatment of the Ligon-Schaaf regularization map, but their calculations 
are still more laborious than one would like [3J. 
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The main novelty of this article is to show that the Ligon-Schaaf regu- 
larization map is almost trivially understood as an adaptation of the Moser 
regularization map. For this reason we recall in Section 2 the Moser regular- 
ization, essentially following the original paper by Moser. The Ligon-Schaaf 
map is the natural adaptation of the Moser map intertwining the Kepler flow 
on the negative energy part P_ of the phase space M? n and the (geodesic) 
Delaunay flow on the punctured cotangent bundle T x of the sphere § n in a 
canonical way. 

The Kepler problem has hidden symmetry in the sense that the visible 
symmetry by the orthogonal group of size n extends to size (n + 1 ) . Of course 
this hidden symmetry becomes naturally visible on § n , and all one has to 
do is to calculate the pull-back under the Ligon-Schaaf map of all angular 
momenta on T x . So our discussion in Section 4 of this hidden symmetry 
is based on the (defined) equivariance of the Ligon-Schaaf map for the full 
symmetry group. In this way one is naturally led to the components of the 
Lenz vector (divided by the square root of —2H) for the extension of Lie 
algebras from so(n) to so(n + 1). 

All arguments of this paper have natural adaptations from the negative 
energy part P_ of phase space to the positive energy part P + of phase space, in 
which case T x becomes the punctured cotangent bundle of hyperbolic space 
HP. Details are easily filled in and left to the reader. At another occasion 
we hope to discuss the implications of the Ligon-Schaaf regularization in 
classical mechanics for the quantum mechanics of the hydrogen atom [12], 

EMU- 

2 Regularization after Moser 

Let us describe the transformation formulas for stereographic projection. Let 
n = (0, • • • , 0, 1) be the north pole of the unit sphere S n in M n+1 , and let 
u = (ui, • • • ,u n+ \) be another point of E> n . Let v = (vi, ■ • • ,v n+ i) be a 
covector on S n at u, so that 



i j 



u^n, u ■ u = 1 , u ■ v = 

are the constraints on u, v in M n+1 . The line through n and u intersects the 
hyperplane MJ 1 orthogonal to n in the point x, and let y be a covector on K n 
at x. The map x i— y u is called stereographic projection, and u i— > x inverse 
stereographic projection. 




Theorem 2.1. The transformation formulas for the inverse stereographic 
projection T*E> n — ¥ T*M n , (u, v) i— >• (x, y) are given by 

Xk = u k /(l - u n+ i) 

Uk = v k (l - u n+1 ) + v n+1 u k 

and for the stereographic projection T*M. n — > T*S n , (x, y) i— > (u, v) become 

u n+ i = (x 2 - l)/(x 2 + 1) 



u k = 2x k /(x 2 + 1) 

Vk = (x 2 + l)y k /2 - (x • y)x k 



v n+1 = x • y 



for k = 1, • • • ,n. These transformations are canonical in the sense that the 
symplectic forms Y^\ dxkAdyk and (the restriction of) J^™ dukAdvk match. 

Proof. Clearly the last coordinate of x = Au + (1 — A)n vanishes, and so 
A(l — u n+ i) = 1. Conversely u = /ix + (1 — fi)n has unit length, and so 
fj,(x 2 + 1) = 2. Hence the transformation formulas between S n and M™ are 
obvious. 

By rotational symmetry we may suppose that n—1. In that case 

xi = «x/(l - u 2 ) , yi = «i(l - u 2 ) + v 2 u x 



and the canonical one form on T*M becomes 

yidxi = (vi(l - u 2 ) + v 2 u 1 )(du 1 /(l - u 2 ) + u 1 du 2 /(l - u 2 ) 2 ) 

= vidui + ViUidu 2 /(l - u 2 ) + v 2 Uidui/(l - u 2 ) + v 2 u\du 2 /(l - u 2 ) 2 
= vidui + v 2 {-u 2 du 2 + uidui + (1 + u 2 )du 2 }/(l - u 2 ) 
= Vidui + v 2 du 2 + v^uxdux + u 2 du 2 }/(l - u 2 ) 
= vidui + v 2 du 2 

as should, using viUi + v 2 u 2 = 0, -uf + u\ = 1, u\dui + u 2 du 2 = 0. 
For the inverse mapping again suppose n = 1. Then we have 

2/1^1 = W 1 - u 2) + ^2«l}wi/(l - «2) = ^1«1 + ^2(1 + U 2 ) = V 2 

using u"l + ul = I, V1U1 + v 2 u 2 = 0, and therefore 

vi = (yi - v 2 U!)/(l - u 2 ) = yi(x 2 + l)/2 - (yiXi)xi 

as desired, using (1 — u 2 ) = 2/(x 2 + 1),-U!/(1 — u 2 ) = x\. Since any diffeo- 
morphism of manifolds induces a symplectomorphism of cotangent bundles 
the theorem follows. □ 

Corollary 2.2. The invariant metric on S n transforms to 

v * = (x 2 + 1) V/4 

under stereographic projection. 
Proof. Indeed the relation 

Vk = {x 2 + l)yk/2 - (x • y)x k , v n+1 = x • y 
for k — 1 • • • , n gives 

v 2 = {x 2 + l)V/4 + (x • y) V - (x 2 + l)(x • y) 2 + (x • y) 2 

and the formula follows. □ 

Consider the function -F(u, v) = v 2 /2 on the cotangent bundle T*S n of 
the sphere § n . The trajectories of the Hamiltonian flow of F on the level 
hypersurface F — 1/2 in T*S n project onto great circles on S n traversed in 



arc length time s with period 2tt. By the above corollary the trajectories of 
the Hamiltonian flow on T*M n of the function 

F(x,y) = (x 2 + l)V/8 

are images under stereographic projection of the trajectories of this geodesic 
flow on T*§ n . Consider the two related functions 

G(x, y) = y/2F(x,y) - 1 = (x 2 + l)y/2 - 1 
H (x, y) = G(x, y)/y - 1/2 = x 2 /2 - l/y 

on T*~R n . The Hamiltonian vector fields of F and G on T^IR™ coincide on the 
level hypersurface F — 1/2 (or equivalently G = or iT = —1/2), since the 
derivative 

{r h^ v^ - 1}' = {r H- 1/^27} 

is equal to 1 for r = 1/2. Hence we have proved the following regularization 
theorem of Moser [10] . 



Theorem 2.3. On the level hypersurface F = 1/2 the trajectories of the 
Hamiltonian flow of the function F(u,\) = v 2 /2 on T*(§™ — {n}) traversed 
in time s equal to arc length transform under inverse stereographic projection 
to trajectories of the Hamiltonian flow of the function H (x, y) = x 2 /2 — l/y 
traversed in real time t on the level hypersurface H = —1/2 with 

ds 1 
dt y 

relating the two time parameters s and t. 

If we set x = p, y = — q then (q, p) \-¥ (x, y) is a canonical transformation 
of R 2n , called the geometric Fourier transform, and 

H(q,p) = H(x,y)=p 2 /2-l/q 

becomes the Kepler Hamiltonian. 

Corollary 2.4. The Moser regularization map $m : T*M. n — >■ T*(§ n — {n}) is 
defined as the composition of stereographic projection with geometric Fourier 
transform. It is a symplectomorphism and explicitly given by the formula 
(q,p) ^ $ M (q,p) = (u,v) with 

u = (2p/(p 2 + 1), 2p 2 /(p 2 + 1) - 1) , v = (-(p 2 + l)q/2 + (q • p)p, -q ■ p) 



as is clear from Theorem \2.1[ On the level hypersurface H = — 1/2 the Moser 
map transforms the Kepler flow with time t to the geodesic flow on E> n — {n} 
with arc length time s, with the infinitesimal Kepler equation 

ds 1 
dt q 

relating the two different time parameters. 

In this corollary the equivalence of the Kepler problem with the geodesic 
flow of the sphere is established only on the energy hypersurface H = —1/2. 
The general case of negative energy H < can be reduced to this case by a 
scaling of variables with 1R + = {p > 0} according to 

q H- p 2 q , p H- p _1 p , H H- p~ 2 H , t H- p 3 t 

and the symplectic form u = J^ dq^ A dp^ scales according to w ^ pu. This 
ends our discussion of the Moser regularization. 

3 Regularization after Ligon and Schaaf 

Let us fix some notations. The phase space T*M n of the Kepler problem with 
canonical coordinates (q, p) will be denoted by P. Let us denote 

P- = {(q 1 p)eP;q^o,%p)<0} 
P-l/2 = {(q,p)eP;q^o,tf(q,p) = -l/2} 

with H (q, p) = p 2 /2 — \jq the Kepler Hamiltonian as before. Likewise let 
the cotangent bundle T*§ n = {(u,v) G T*R n+1 ;u = l,u ■ v = 0} of the 
sphere will be denoted by T, and put 

T x = {(u,v)eT;v^0} 

T_ = {(u,v) eT;u^n,v^0} 

T-i/2 = {(u,v)GT]u/n,t) = l} 

for the various submanifolds of T. 

Definition 3.1. The Delaunay Hamiltonian H is defined by 

H(u,v) = ~ 

2v z 

on the punctured cotangent bundle T x ofS n . 
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It is clear that the flow of the Delaunay Hamiltonian has a similar scaling 

symmetry 

uh->u,vH-pv,ffH' P~ 2 H, t (->■ p 3 t 

as the Kepler problem, and likewise the symplectic form u = J2 duu A dvk 
scales according to u) h-> pa). 

Remark 3.2. 27ie Moser fibration Um '■ P — > P-1/2 with Um (q, p) = (u, v) 

is defined by 

u = (y/-2Hqp, qp 2 - 1) , v = (-q/g + (q • p)p, -V~2Hq • p) . 

iVoie £/iat Iljv/ is invariant under the action of the scale group R + on P_ . On 
the submanifold P-1/2 we have (p 2 + l)/2 = 1/q, and therefore 

$A/(q,p) = n M (q,p) 
on P_i/2 fry a direct comparison with the formula for the Moser map $^ 



m Corollaru\2.4\ In turn this equality implies that Um is indeed a smooth 



fibration with fibers the orbits of the scale group R + . 

Definition 3.3. The Ligon-Schaaf regularization map §ls : P- — > T_ zs 

defined by $1,5 (q, p) = (r, s) itn't/i 



r = ((cosw n+ i)u+ (sinw n+ i)v) , s = ((- sinw n+ i)u + (cosw n+ i)v)/v / -2P 
with u,v6§"C M n+1 gzven 6?/ 

u = (v /= 2# g p, gp 2 - 1) , v = (-q/g + (q • p)p, -^^2H(q • p)) 
as the components of the Moser fibration IlM(q, p) = (u, v). 

By definition the Ligon-Schaaf map is equivariant for the action of the 
scale group 1R + on P_ and T x , which in turn is equivalent to $2s-^ = H. 

For (u, v) G XI 1/2 the vectors u, v form an orthonormal basis in the plane 
Ku + Kv. If we define a complex structure on this plane by 

IVL = V , iv = — u 

then the Ligon-Schaaf map can be written in the compact form 



r = e lVn+x u , s = e Wn+1 v/V-2H . 



We claim that for (q, p) G P-1/2 and (u, v), (r, s) G IL1/2 related as above 
the following three differential equations 

rfq dp 3 

du dv 

ds ds 

dv ds 

dt dt 

are equivalent with time parameters t and s related by ds/dt = 1/q. The 
first two of these are equivalent under the Moser map $m by Corollary 12.41 
The equivalence with the third one follows from 

dr dv n+ \ _•„ dsdu , 9 

* = ! if r + e mTs = <- p + 1/ « )s + s >o = t- 2H > = s 

and likewise for s = ir. This shows that the first and the third equation 
are equivalent under the Ligon-Schaaf map §ls- Using that the Kepler and 
the Delaunay Hamiltonians have the same scale symmetry and that &ls is 
equivariant for the two scale symmetries it follows that $^5 intertwines the 
Kepler flow on P_ and the Delaunay flow on T_. 

Since $^5-^ = H it follows that $^5 : P_ — > T__ is a symplectomorphism. 
Indeed $>ls is obtained from the canonical Moser map $m : P — > T by a 
canonical modification along the ruled surfaces with base the Kepler orbits 
in P-1/2 and rulings the scale action. For such a ruled surface S in P_ the 
form u\s = dt A dH is mapped under $^5 to the form £j\g = dt A dH on 
S = ^ls(S). This proves the following result of Ligon and Schaaf [8j. 

Theorem 3.4. The Ligon-Schaaf regularization map $^5 : P_ — > T_ is a 
symplectomorphism intertwining the Kepler flow on P_ with the Delaunay 
flow on T_ (with respect to the same time parameter t). 

The regularization is obtained by the partial compactification T_ c — )■ T x . 
The incomplete Kepler flow of the collision orbits on P_ is regularized under 
<£>£S by T_ ■— y T x to the complete Delaunay flow on T x . For this reason the 
punctured cotangent bundle T x of § n is called the Kepler manifold [13]. We 
now move on to discuss the symmetry of the Kepler problem. 



4 The symmetry group 

If we write Ly = qiPj — q^pi for the components of angular momentum 
L = q A p on the phase space R 2n then {Lij, Ljk} = L^i while {Ly, Lj-i} = 
if #{i,j, k, 1} j^ 3. These are the commutation relations for the Lie algebra 
so(n). 

Theorem 4.1. On the negative energy part P_ of phase space put 

Li( n +i) = —L(n + x)i = K i /\ / —2H 

with Ki the components of the Lenz vector 

K = (p 2 - l/g)q- (q-p)p. 

If we write M^ = riSj — rjSi for the components of the angular momentum 
M = r A s on the phase space R 2 ( n+1 ) then 

for i,j = 1, ■ •• ,n + 1. In particular {Lij,Ljk} = Lki and {Lij,Lki} = 
if #{i,j,k,l} 7^ 3 for i,j,k,l = 1, ••• ,n + l. The conclusion is that the 
Ligon-Schaaf regularization map 

intertwines the two infinitesimal Hamiltonian actions of so(n + 1). 

Proof. Using the formula for $^5 of Definition 13. 31 the relation <§>* LS Mij = L^ 
is a rather straightforward exercise in substition. 

Indeed for i,j = 1, • • • , n and My = TjSj — TjSj the expression for 3>* LS Mij 
a t (q, p) is equal to 



{(cosv n+1 )(V-2Hqpi) + (sin v n+1 )(-qi/q+ (q-p)p;)} 
x{(-smv n+1 )(\ / -2Hqp j ) + (cos v n+1 )(-qj/q + (q • p)pj)}/y/-2H 

{(cosv n+ i sinv n+1 )(-q 2 p i p j )}V-2H 

+{(cos 2 v n+1 )(-piqj + q(q ■ v)PiPj) + (sin 2 v n+1 )(qiPj - g(q ■ v)PiPj)} 

+{(cosw n+ i smv n+1 )(-qi/q + (q • p)pi){-qj/q + (q • p)pj)}/y/-2H 

-{* ++ j) 

{cos 2 v n+1 (qiPj - q^) + sin 2 v n+1 (qiPj - qjPi)} = (qiPj - qjPi) = L tj 



with {i <H> j} denoting the same expression with % and j interchanged. Of 
course this formula is also obvious because our geometric construction of the 
Ligon-Schaaf regularization map is equivariant for the action of so(rz). 

Likewise for % — 1, • • • , n and Mj( n+1 ) = TjS n+ i — r n+ iSi the expression for 
®ls M i(n+i) at (q, p) is equal to 



= {(cosv n+1 )(\ / -2Hqp i ) + (sin v n+1 )(-qi/q + (q • p)Pi)} 

x{(- smv n+1 )(qp 2 - 1) + (cosv n+1 )(-V-2H(q ■ p))}/V^2H 
-{(cosv n+1 )(qp 2 - 1) + (smv n+1 )(-y/-2H(q • p))} 
x{(- sin v n+1 )(V-2Hqpi) + (cos v n+1 )(-qi/q + (q- p)pi)}/y/-2H 
= {(-cosf n+ isinf n+ i + cosv n+ isinv n+1 )(qp 2 - l)qpi] 

+{(- sinw n+1 cosw n+ i + sinw n+ i cosv n+ i)(q • p)(-ft/g + (q • p)Pi)} 
+{(cos 2 w n+ i + sm 2 v n+ i)(2Hq(q-p)pi)}/V-2H 
+{(-sin 2 w n+ i -cos 2 v n+1 )(qp 2 - l)(-ft/g + (q • p)p i )}/\ / -2H 
= {2Hq(c L ■ p) Pi - (qp 2 - l)(-q t /q + (q • p) Pi )}/V=2H 
= {(QP 2 " 2)(q • p) Pl + (qp 2 - l)qi/q - (qp 2 - l)(q • p^j/v 73 ^ 
= {(p 2 - l/g) ft - (q • p) Vl }/^2H = K t /^2H = L i{n+1) 

which proves the theorem. D 

Corollary 4.2. On P_ the Hamiltonian vector fields of the functions L^ 
for i,j = 1, • • • ,n + l integrate to an incomplete Hamiltonian action of the 
rotation group SO(n + 1). Moreover the Ligon-Schaaf map 

intertwines this action with the standard complete Hamiltonian action of 
SO(n + 1) on the punctured cotangent bundle T x of'E> n . In particular $ls 
intertwines the corresponding moment maps 

\x: P -»■ so(n + 1)* , \x : T x ->■ so(n + 1)* 

in £/ie sense £/ia£ /i(q, p) = //(<&£g(q, p)). In £nrn £/ws implies that 

fi 2 = L 2 + K 2 /(-2H) = l/(-2H) 

because ft 2 = l/(—2H). 
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The incomplete Hamiltonian action of so(n+l) on P_ is regularized under 
&ls by the partial compactification T_ ■=->■ T x to a complete Hamiltonian 
action of so(n + 1) on T x . The fact that &ls intertwines the two moment 
maps was the point of departure of Cushman and Duistermaat [3]. The next 
corollary is due to Gyorgyi [5] , but on the more refined quantum level it was 
obtained before by Pauli in his spectral analysis of the hydrogen atom 



Corollary 4.3. On the phase space {(q, p) G P; q ^ o} of W l — {o} the 

components Ki of the Lenz vector satisfy the Poisson bracket relations 

{Lij, K k } = b ik Kj - 5 jk Ki , {K^ Kj} = -2HL i;j 

fori,j,k = l,-- ,n. 

Proof. These Poisson brackets hold on P by the above theorem, and extend 
from P to {(q, p) 6 P; q ^ o} as analytic identities. □ 

In turn this implies that on the positive energy part P + of phase space 
the symmetry algebra becomes the Lorentz algebra so(n, 1) rather than the 
orthogonal algebra so(n + 1). 

References 

[1] V. Bargmann: Zur Theorie des Wasserstoffatoms, Bemerkungen zu der 
gleichnamigen Arbeit von V. Fock, Zeitschrift fur Physik 99 (1936), 576- 
582. 

[2] Bruno Cordani: The Kepler Problem, Progress in Mathematical Physics 
29, Birkhauser, Basel, 2003. 

[3] R.H. Cushman and J.J. Duistermaat: A Characterization of the Ligon- 
Schaaf Regularization Map, Comm. Pure and Appl. Math. 50 (1997), 

773-787. 

[4] Victor Guillemin and Shlomo Sternberg: Variations on a Theme by 
Kepler, AMS Colloquium Publications, Volume 42, 1990. 

[5] G. Gyorgyi: Kepler's equation, Fock variables, Bacry's generators and 
Dirac brackets, Nuovo Cimento 53A (1968), 717-735. 



11 



[6] 
[7] 
[8] 
[9] 
[10 

[11 



[12 



[13 



V. Fock: Zur Theorie des Wasserstoffatoms, Zeitschrift fur Physik 98 
(1935), 145-154. 

Maris van Haandel and Gert Heckman: Teaching the Kepler Laws for 
Freshmen, Math. Intelligencer 31:2 (2009), 40-44. 

Thomas Ligon and Manfred Schaaf: On the Global Symmetry of the 
Classical Kepler Problem, Reports on Math. Phys. 9 (1976), 281-300. 

John Milnor: On the Geometry of the Kepler Problem, Amer. Math. 
Monthly 90:6 (1983), 353-365. 

J. Moser: Regularization of Kepler's Problem and the Averaging Method 
on a Manifold, Comm. Pure and Appl. Math. 23 (1970), 609-636. 

I. Newton: The Principia, Mathematical Principles of Natural Philoso- 
phy, A New Translation by LB. Cohen and A. Whitman, with A Guide 
to Newton's Principia by LB. Cohen, University of California Press, 
Berkeley, 1999. 

W. Pauli: Uber das Wasserstoffspektrum vom Standpunkt der neuen 
Quantenmechanik, Zeitschrift fur Physik 36 (1926), 336-363. 

J.-M. Souriau: Sur la variete de Kepler, Symp. Math. 14 (1974), 343-360. 



Gert Heckman, Radboud University Nijmegen, P.O. Box 9010, 
6500 GL Nijmegen, The Netherlands (E-mail: g.heckman@math.ru.nl) 
Tim de Laat, Radboud University Nijmegen, P.O. Box 9010, 
6500 GL Nijmegen, The Netherlands (E-mail: tlaat@science.ru.nl) 



12 



